A contractor-renormalization study of Hubbard plaquette clusters 



O 
(N 

> 

O 






I 

O 

o 



> 

oo 
l> 

p 

o 
o 



X 



Shirit Baruch and Dror Orgad 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

(Dated: November 4, 2010) 

We implement the contractor-renormalization method to study the checkerboard Hubbard model 
on various finite-size clusters as function of the inter-plaquette hopping t' and the on-site repulsion 
U at low hole doping. We find that the pair-binding energy and the spin gap exhibit a pronounced 
maximum at intermediate values of t' and U, thus indicating that moderate inhomogeneity of the 
type considered here substantially enhances the formation of hole pairs. The rise of the pair- 
binding energy for t' < imax is kinetic-energy driven and reflects the strong resonating valence bond 
correlations in the ground state that facilitate the motion of bound pairs as compared to single 
holes. Conversely, as t' is increased beyond i^ax antiferromagnetic magnons proliferate and reduce 
the potential energy of unpaired holes and with it the pairing strength. For the periodic clusters 
that we study the estimated phase ordering temperature at t' = t^ax is a factor of 2-6 smaller than 
the pairing temperature. 

PACS numbers: 74.81.-g, 71.10.Fd 



I. INTRODUCTION 

It is by now generally accepted that spatial inho- 
mogeneity may emerge either as a static or as a fluc- 
tuating effect in strongly-coupled models of the high- 
temperature superconductors, and indeed in many of the 
real materials ii What is far from being settled is the issue 
of whether such inhomogeneity is essential to the mecha- 
nism of high-temperature superconductivity from repul- 
sive interactions. While most researchers would proba- 
bly answer this question in the negative one should bare 
in mind the absence of a conclusive evidence that the 
single-band two-dimensional Hubbard model, widely be- 
lieved to be the "standard model" of high-temperature 
superconductivity, actually supports superconductivity 
with a high transition temperature.^ On the other hand, 
when examined on small clusters the same model and 
its strong-coupling descendent, the t ~ J model, exhibit 
robust signs of incipient superconductivity in the form 
of a spin-gap and pair binding.^ This fact points to the 
possibility that the strong susceptibility towards pairing 
is a consequence of the confining geometry itself. 

This line of thought has been pursued in the past by 
considering the extreme limit where the electronic den- 
sity modulation is so strong that the system consists of 
weakly coupled Hubbard ladders^ or plaquettesS. Be- 
yond the questionable applicability of such models to 
the physical systems, which are at most only moderately 
modulated, it is clear that strong inhomogeneity, even if 
beneficial to pairing, is detrimental to the establishment 
of phase coherence and consequently to superconductiv- 
ity. On both counts it is, therefore, desirable to extend 
the analysis to the regime of intermediate inhomogeneity. 

Recently, the checkerboard Hubbard model, con- 
structed from 4-site plaquettes with nearest-neighbor 
hopping t and on-site repulsion [/, was studied as func- 
tion of the inter-plaquette hopping t' (see Fig. [T]). Tsai 
et al^ diagonalized exactly the 4x4 site cluster (2x2 



plaquettes) and found that the pair-binding energy, as 
defined by Eq. ([2]) below, exhibits a substantial maxi- 
mum at t' w i/2 for [/ w 8t and low hole concentration. 
Doluweera et al^, on the other hand, used the dynam- 
ical cluster approximation in the range 0.8 < t' /t < 1 
and obtained a monotonic increase in both the strength 
of the d-wave pairing interaction and the superconduct- 
ing transition temperature, Tc, towards a maximum that 
occurs in the homogeneous model. 

In this paper, we use the contractor-renormalization 
(CORE) method^ to derive an effective low-energy 
Hamiltonian for the checkerboard Hubbard model, which 
we then diagonalize numerically on various finite-size 
clusters. We begin by establishing the region of appli- 
cability of the CORE approximation by contrasting its 
predictions with the exact results of Ref. for 2x2 pla- 
quettes. Our findings indicate that at low concentrations 
of doped holes the two approaches agree reasonably well 
unless t' is larger than a value, which increases with U. 
Deviations also appear for small t' when U is large. We 
identify probable sources of these discrepancies. 

Based on the lessons gained from the small system 
we go on to study larger clusters of up to 10 plaque- 
ttes. These include the periodic 6x6 sites cluster and 
2-leg and 4-leg ladders with periodic boundary condi- 
tions along their length. Within the region where CORE 
is expected to provide reliable results the pair-binding 
energy continues to exhibit a non-monotonic behavior 
with a pronounced maximum at intermediate values of 
t' and U. The precise location of the maximum depends 
on the cluster geometry but it typically occurs in the 
range i^^x ~ 0.5 — 0.7i and [/max ~ 5 — St. The spin 
gap of the doped system follows a similar trend, often 
reaching the maximum slightly before the pair-binding 
energy. These findings demonstrate that moderate inho- 
mogeneity, of the type considered here, can substantially 
enhance the binding of holes into pairs. 

In an effort to elucidate the source of the maximum 



we have looked into the content of the ground state and 
calculated the contributions of various couplings in the 
effective Hamiltonian to its energy. Our results indicate 
that for t' < t'^g,^ the doped holes move in a background, 
which is composed predominantly of plaquettes that are 
in their half-filled ground state. This background pos- 
sesses strong intra-plaquette singlet resonating valence 
bond (RVB) correlations, which facilitate the propaga- 
tion of pairs relative to independent holes. The rise in 
the pair-binding energy while t' grows towards t'^^^^ is a 
result of a faster decrease of the pair kinetic energy in 
comparison to that of unpaired fermions. As t' crosses 
tj-fjjjx s-i^d approaches the uniform limit the ground state 
contains a growing number of plaquettes that support 
antiferromagnetic (AFM) magnons. In this regime of in- 
creasing AFM correlations the kinetic energy changes rel- 
atively little with i', and the decrease of the pair-binding 
energy for t' > ij^^^^ is caused by the lowering of the 
energy of single holes due to their interactions with the 
magnons. Interestingly, we find that the maximum in the 
pair-binding energy of the periodic clusters is accompa- 
nied by a change in the crystal momentum of the single- 
hole ground state from the F — M and symmetry related 
directions at t' < ij^^^ to the Brillouin-zone diagonals at 
t' > ijjjax- ^ similar correlation was also found for the 
3-hole ground state of the 6x6 sites cluster. 

While the pair-binding energy sets a pairing scale, Tp, 
a phase-ordering scale, Tg, is provided by the phase stiff- 
ness. The latter was evaluated from the second derivative 
of the ground state energy with respect to a phase twist 
introduced by threading the system with an Aharonov- 
Bohm flux. We have found that as the twist is taken to 
zero, the CORE energy curvature typically converges to- 
wards a limiting value only when t' < t^^g^^. Within this 
region the phase stiffness increases monotonically with t' . 
Our results indicate that for the lightly doped periodic 
clusters that we have considered phase fluctuations dom- 
inate over pairing, specifically, Tp « 2 — GTg at t' = tj„ax- 
The limitations of the present study make it difficult to 
draw conclusions regarding the behavior of Tc in the two- 
dimensional thermodynamic limit. 

We have also calculated the pair-field correlations be- 
tween Cooper-pairs that reside on the most distant bonds 
allowed by our finite clusters. As expected, these corre- 
lations are consistent with d-wave pairing. However, in 
contrast to the pair-binding energy and the phase stiff- 
ness the correlations change little with t' and are small in 
magnitude. This discrepancy might be resolved in light 
of our finding that only few holes are tightly bound into 
pairs that reside within a single plaquette. Moreover, we 
obtain that the number of such pairs changes relatively 
little with f with no apparent correlation to the substan- 
tial maximum in the pair-binding energy. Taken together 
these findings suggest that the correlation function which 
we and others often use to identify and quantify pair- 
ing in the Hubbard model may be ill-constructed to take 
account of the more extended and structured nature of 
pairing in this model. 
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FIG. 1: The checkerboard Hubbard model. Shown here are 
two of the clusters that we studied. The bonds labeled ab, 
cd, and ef specify locations used in calculating the pairing 
correlations. 



II. MODEL AND METHOD 

The Hamiltonian of the checkerboard Hubbard model, 
which wc have studied, is given by 

^ = ~ E (*y4-Cj.- + H.C.) +UY,n^,t^^,h (1) 



{i,j),<^ 



where cj ^ creates an electron with spin a =t) i s-t site i 

of a two-dimensional square lattice. Here rti^ = cj g.Cj,cr, 
and (i,j) denotes nearest-neighbor sites. The hopping 
amplitude is Uj = t for i and j on the same plaquette, 
while tij = t' when they belong to neighboring plaque- 
ttes, as shown in Fig. [TJ 

The first step in obtaining the CORE effective Hamil- 
tonian for the above model, is the exact diagonalization 
of a four-site plaquette. Out of the full spectrum, the M 
lowest-energy states are retained. The reduced Hilbert 
space, in which the effective Hamiltonian operates, is 
spanned by the tensor products of these states on differ- 
ent plaquettes. Next, the Hamiltonian ([1]) is diagonalized 
on N connected plaquettes and the A/^ lowest-energy 
states are projected onto the reduced Hilbert space and 
Gram-Schmidt orthonormalized. Finally, after replac- 
ing the exact eigenstates by their projections, the A^- 
plaquette Hamiltonian can be represented as one for M 
types of hard core particles coupled via A^-body inter- 
actions. The CORE approximation consists of applying 
the resulting effective Hamiltonian to the study of larger 
clusters. By construction, the spectrum of the CORE 
Hamiltonian coincides with the low-energy spectrum of 
the exact problem on A^ plaquettes. We note, however, 
that this ceases to be the case if one or more of the ex- 
act low-energy states have zero projection on the reduced 
Hilbert space, or, if some of them are projected onto the 
same tensor-product state. In the following we demon- 
strate that such a problem arises in certain parameter 
regions of the model dH). 

We concentrate on relatively low hole densities as mea- 
sured from the half-filled system. The simplest trunca- 



tion used to describe this regime is to retain the ground 
state of the half- filled plaquette [a total spin singlet 5 = 
with plaquette momentum q — (0,0)], its 5 = f, q = 
(tt, tt) triplet of lowest lying AFM magnon excitations, 
and the S" = 0, q = (0, 0) hole pair ground state. ^ The in- 
clusion of the magnon excitations is essential for retriev- 
ing the correct magnetic behavior at low hole dopingiiS 
Below we show that they also play an important role in 
the physics of hole binding. One can improve the approxi- 
mation by including in the CORE plaquette basis also the 
two degenerate doublets S^ — ±1/2, q = (0, tt), (tt, 0), 
comprising the single hole ground state .-'-'-'i"'"^ Moreover, 
the inclusion of these states is mandatory for the purpose 
of calculating the pair binding energy, which is one of 
the goals of the present work. Consequently, our CORE 
scheme consists of keeping the above mentioned M = 9 
states. We have considered only range-2 interactions, i.e. 
TV = 2. 

The resulting effective Hamiltonian includes all pos- 
sible couplings, which respect the symmetries of the 2- 
plaquettes problem. These include the conservation of 
number of holes Nh, invariance under SU(2) spin rota- 
tions and under reflections about the central bonds of 
the cluster in the x and y directions. The latter, to- 
gether with the conservation of N^, imply that within 
our reduced Hilbert space, as defined above, the total 
plaquette momentum qi 4- q2 is also conserved (modulo 
27r). We will not list here the 45 couplings which are al- 
lowed by the symmetries. Instead, we will describe the 
most important ones in the appropriate context and refer 
the reader to the Appendix for a detailed description of 
the Hamiltonian. 

Many of the results reported in the following are de- 
rived from the spectrum of the effective Hamiltonian as 
obtained by exact diagonalization. We also calculate var- 
ious ground-state correlations. To this end we project the 
appropriate operators on the reduced Hilbert space^ be- 
fore evaluating their ground-state correlation function. 



III. RESULTS 



Although the size of the Hilbert space is massively re- 
duced by the CORE approximation it still grows expo- 
nentially with the size of the system. Therefore, even 
the largest clusters that we are able to diagonalize us- 
ing this method are too small for a direct calculation of 
Tc. Instead we calculate various properties of the system 
which are indicative of the two necessary ingredients for 
superconductivity: pairing and phase stiffness. We begin 
with the former and study its behavior as function of t' 
and U on various geometries. These include the 4x4 
and 6x6 periodic clusters, seen in Fig. [I] as well as 2- 
leg and 4-leg ladders with periodic boundary conditions 
along their length, which extends up to 20 sites. 




FIG. 2: The pair-binding energy in a periodic 4x4 cluster 
at 1/16 hole doping as obtained by (a) CORE, and (b) ex- 
act diagonalization (Ref. 6). CORE projects out low energy 
states from the effective Hilbert space in the region above the 
dashed line. The crystal momentum of the degenerate single- 
hole ground state is (0,7r) and (7r,0) below the solid line and 
(0,0) and (7r,7r) above it. 




FIG. 3: The pair-binding energy in a periodic 4x4 cluster at 
1/16 doping for various values of the interaction strength (a) 
U = 4f, (b) U = 8t, and (c) U = lOt. Triangles depict the 
CORE results and circles correspond to the exact diagonal- 
ization results of Ref. [3. 



A. Pair-binding energy and spin-gap 

The pair-binding energy is defined by 

Apb{M/N) = 2Eo{M)-[EoiM + 1) + EoiM - 1)] , (2) 



where Eq{M) is the ground-state energy of the system 
with M holes doped into the A^-site half-filled cluster. 
Consider two identical clusters each with M holes. If 
holes tend to pair and M is odd it should be energetically 
favorable to move an electron from one cluster to another 
in order to obtain a fully-paired state in both. On the 
other hand, such a redistribution should be unfavorable 
if M is even. In this sense, a positive Apb for odd M and 
a negative Apj, for even M signifies an effective attraction 
between holes. 

Recently, Tsai et al^ have found by exact diagonal- 
ization of the periodic 4x4 cluster that the pair-binding 
energy exhibits a pronounced maximum both as function 
of t' and U . Their results allow for a critical evaluation 
of the validity of the CORE method in a range of pa- 
rameters. To this end we present in Figs. [2] and [3] a 
comparison between the CORE and the exact results for 
Apb(l/16). It is clear that CORE introduces substan- 
tial errors in two specific regimes: small U and large t' 
[Fig. El^a)], and large U and small t' [Fig. El^c)], while it 
is in reasonable agreement with the exact results in the 
intermediate parameter regime. 

An obvious source for the discrepancies is the fact 
that our CORE approximation includes only range-2 cou- 
plings. Longer-range interactions are expected to become 
more important as the system becomes more homoge- 
neous i.e. when t' —?' t. We believe that the devia- 
tions between the CORE predictions and the exact re- 
sults in this limit, especially for small U where the pair 
size is expected to be large, are mainly due to insuffi- 
cient range of the effective interactions. A related prob- 
lem may emerge at large U where the extent of magnetic 
correlations grow. However, we did not confirm these 
conjectures by explicit calculations. 

A more subtle source of errors, which we have men- 
tioned already in the previous Section, is the fact that 
low-energy states may be projected out from the CORE 
effective Hilbert space in the process of generating the 
effective Hamiltonian. This happens when a low-lying 
state of a connected cluster has zero overlap with the 
tensor-product states of the effective Hilbert space or 
when two or more low-lying states are mapped onto the 
same state in the effective space (Note, however, that 
spin-rotation symmetry is preserved in the sense that 
spin multiplets are either kept or projected out as a 
whole.) Fig. U] depicts for each of the sectors in which 
such a problem arises the excitation energy of the lowest 
projected-out state in units of the bandwidth of the kept 
states in the sector. We also denoted in Figs. [5] and [5] 
the parameter region where the problem occurs. 

The overlap issue is responsible for the failure of CORE 
in the regime of small t' and large U. When U > 7.858i 
and for t' — the N^ — I, S = 3/2 double-plaquette 
(eight- fold degenerate) ground state \Nh = 1,5 = 3/2)2 
consists of one plaquette in its half- filled ground state and 
a second plaquette in a fully-polarized S = 3/2 single- 
hole state. The latter resides outside the effective Hilbert 
space and therefore \Nh = 1,5' = 3/2)2 is projected out. 
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FIG. 4: The excitation energy of the lowest energy state that 
is projected-out by CORE in units of the bandwidth of the 
kept states in its sector. 



This ceases to be the case once t' is turned on as a result 
of a component which appears in \Nh = 1,S — 3/2)2 
and corresponds to a system with a magnon on one pla- 
quette and a plaquette-fermion on the other. However, 
the amplitude of this component diminishes with increas- 
ing U. This leads CORE to misidentify the nature of 
\Nh = 1,5 = 3/2)2 and induces an abrupt increase in 
the magnon-fermion interaction [Vf^ '"'*' in Eq. (|A5|) ] 
for small t' . As a result, CORE underestimates the en- 
ergy of the two-hole ground state of the 4x4 cluster and 
consequently predicts an erroneously large pair-binding 
energy, see Fig. [31 Nevertheless, it appears that away 
from this region of parameters the projected-out states 
are high enough in energy as to not cause qualitative 
errors. 



Based on the comparison of Apt depicted in Figs 
and similar plots presented below for the spin-gap [Fig. 
El^a)] and pair-field correlations [Fig. [T^ a)] we conclude 
that CORE agrees semi-quantitatively with the exact re- 
sults provided U/50 < i' < U/8. Within this region, and 
across all geometries studied, we found the pair-binding 
energy to exhibit the same qualitative behavior consist- 
ing of a broad peak both as function of t' and U. This 
conclusion holds true also when one varies the doping 
level (at least in the low-doping regime which we have 
considered) as can be seen from the results for the 6x6 
cluster presented in Fig. [S] In addition, the same fig- 
ure suggests that the above mentioned problems with 
the CORE method become less severe as the size of the 
system increases. 

The association of positive pair-binding energy with 




FIG. 5: The pair-binding energy in a periodic 6x6 cluster 
at (a) 1/36, and (b) 3/36 hole doping. CORE projects out 
low energy states from the effective Hilbert space in the region 
above the dashed line. In (a) the crystal momentum of the de- 
generate single-hole ground state is (0, ±27r/3) and (±27r/3, 0) 
below the solid line and (±27r/3, ±27r/3) above it. In (b) the 
crystal momentum of the degenerate 3-hole ground state is 
(±27r/3, ±27r/3) between the solid lines and (0, ±27r/3) and 
(±27r/3, 0) elsewhere. 




FIG. 6: The spin gap of undoped and two-hole doped systems 
at [/ = 8i. (a) The 4x4 periodic cluster - comparison between 
CORE and exact diagonalization results, (b) CORE results 
for the spin-gap As and the pair-binding energy Ap(,(l/36) of 
the 6x6 periodic cluster. 



Cooper pairing may be contested on the ground that it 
can also be taken as evidence for a tendency of the sys- 
tem to phase separate. We beheve that this is not the 
case for the model studied here for the following rea- 
sons. First, in accordance with the interpretation dis- 
cussed above of Aph as indication for hole pairing we have 
found its sign to change according to (—1)*^+^ for all the 
clusters and doping levels which we have considered. Sec- 
ond, while the appropriate criteria for identifying regimes 
of phase separation from finite size studies include the 
Maxwell constructior^i^ and measurements of the sur- 




FIG. 7: The spin-gap As and pair-binding energy Apt at 
U = 8t for two-hole doped a) 2-leg ladders, (b) 4-leg ladders. 



face tension in the presence of boundary conditions that 
force phase coexistence, a crude way of identifying phase 
separation is by calculating the inverse compressibility 
K^^ — n^d^/dn, where /i is the chemical potential and n 
the electronic density. For numerical purposes a discrete 
version is used, which in our case reads 

K'^ ex Eo{M + 2)+ Eo{M - 2) - 2Eo{M). (3) 

Negative inverse compressibility indicates instability to- 
wards phase separation. We always find k^^ > 0. Fi- 
nally, whenever the ground state is a spin singlet one 
can define the spin-gap as the energy gap to the lowest 
5 = 1 excitation. We have calculated the spin gap for 
the two-hole doped systems and found that in all cases it 
follows the pair-binding energy in the regime of small to 
moderate t', see Figs. [6]and[7l This coincidence strongly 
suggests that in this regime the lowest S = 1 excitation is 
a result of a dissociation of a hole pair into two separate 
holes. It is interesting to note that we always observe 
that the spin-gap reaches a maximum and starts to drop 
before the pair-binding energy does so. This may be 
an indication that moderate inhomogeneity supports the 
formation of a bound 5 = 1 magnon-hole-pair state J^ii^ 
Consequently, our findings and the above arguments 
lead us to conclude that inhomogeneity of the type in- 
cluded in the checkerboard Hubbard model substan- 
tially enhances hole-pairing. The precise position of the 
point of optimal inhomogeneity in the sense of strongest 
pairing depends on the cluster geometry and interac- 
tion strength. Albeit, it typically occurs in the range 
t' sa 0.5— 0.7t and C/max ~ 5— 8t. We note that this fact 



implies that the physics behind the large pairing scale of 
the model necessarily involves inter-plaquette couplings 
since the single plaquette does not support hole-pairing 
beyond Uc ~ iM^ 



B. Energetics and structure of the ground state 

What drives the enhancement of hole-pairing and what 
is the reason for its maximum as function of t'? In an 
attempt to gain insights into these questions we have 
took advantage of the fact that CORE provides us with 
an effective Hamiltonian whose various couplings can be 
classified and analyzed. To this end we have divided 
the 45 different couplings into four groups, as described 
in the Appendix. They include; fermion and hole-pair 
"bare" kinetic terms (including fermion and pair hop- 
ping as well as Andreev-like pair creation and disin- 
tegration), magnon-assisted fermion and pair hopping, 
fermion and pair interactions and finally, interactions 
involving magnons. Fig. |S] depicts the contribution of 
each group to the ground-state energy of the N^ = 0, 1, 2 
doped 6x6 periodic cluster and to its pair-binding energy 
Ap6(l/36) at U = 8t. 

Fig. [S] makes it clear that the increase in the pair- 
binding energy from t' — to t'^^^^ is dominated by a 
faster decrease of the kinetic energy of hole pairs as com- 
pared to unpaired holes. Furthermore, in this region the 
pair-binding energy is largely determined by the "bare" 
kinetic terms while the (negative) contribution of hop- 
ping processes that involve magnons is much smaller. 
The small contributions of the various interactions ap- 
proximately cancel out. Looking more closely at the way 
charges propagate in this range of t' we found that the 
main channel for single holes is a direct hop between 
neighboring plaquettes but that this process is virtually 
non-existent for hole pairs. Instead, a pair propagates 
predominantly by Andreev-like dissociation into single 
holes on adjacent plaquettes and recombination of these 
holes into a pair one register away from its original posi- 
tion [as described by the last term in Eq. (jA2[) ]. 

For t' > ijjiax the behavior changes qualitatively and 
rather abruptly. The gain in kinetic energy of the pair rel- 
ative to that of unpaired holes ceases to increase. While 
pairs continue to propagate mainly via a series of dissoci- 
ation and recombination events, single holes move almost 
exclusively by hopping processes involving magnons [the 
second and third terms in Eq. (jASP J. The decrease in 
the pair-binding energy in this regime is induced by a 
sharp decrease of the potential energy of the unpaired 
holes owing to their interactions with the magnons. On 
the other hand, the contribution of interactions not in- 
volving the magnons to the pair-binding energy does not 
show a significant change as t' is driven through i^^x- 

The above results suggest that the AFM magnons play 
an important role in inducing the change in the behav- 
ior of the pair-binding energy. To further test this con- 
clusion we have looked at the evolution of the ground- 
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FIG. 8: Ground state expectation values of various effec- 
tive couplings for the 6x6 periodic cluster at U — St: (a) 
fermion and pair hopping; (b) fermion and pair magnon- 
assisted hopping; (c) fermion and pair interactions; (d) in- 
teractions involving magnons; (e) the full Hamiltonian. The 
insets show the contribution of each group of couplings to 
the pair-binding energy. The full binding energy reaches a 
maximum at t' — 0.6t as indicated by the dotted line. 



state content with t' . Fig. IH] shows the average num- 
ber of magnons, fermions and pairs in the Nh = — 4 
ground states of the 6x6 periodic cluster. Evidently, the 
magnons begin to proliferate slightly before the maxi- 
mum in the pair-binding energy is reached. Concomi- 
tantly, there is an increase of AFM correlations in the 
system as can be seen from Fig. [TUl which depicts the 
staggered magnetization m(^ ^r) defined by 



*(;T,7r) 



N 



N ^ J 



(4) 



where Q — (tt, tt) and Sj is the electronic spin operator 
on site j at position rj. In contrast to the behavior of 
the magnons, the fermions-to-pairs ratio does not change 
considerably at moderate values of t'. Note that at i' = 
all the holes appear as single fermions. This is a man- 
ifestation of the absence of pair-binding on the single 
plaquette at U = 8t. 

We have found that the same behavior, both in terms 
of energetics and structure of the ground state, persists 
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FIG. 9: The average number of (a) magnons, (b) fermions, 
and (c) hole pairs in the ground state of the 6x6 periodic 
cluster ed, U — 8t. The position of the pair-binding energy 
maximum is indicated by the dotted line. 



across the entire range of geometries and doping levels 
which we have studied. Therefore, we conclude that the 
initial rise of the pair-binding energy for t' < ij^ax i^ 
kinetic-energy driven. In this range most of the plaque- 
ttes are in their half-filled, RVB-correlated ground-state. 
This type of background facilitates the motion of bound 
pairs as compared to single holes. When t' approaches 
t'^g^ the undoped background changes its nature and be- 
comes more AFM. The gain in kinetic energy associated 
with hole-pairing saturates and instead a gain in the po- 
tential energy of unpaired holes sets in due to their in- 
teractions with the AFM magnons. This leads to the 
decrease of the pair-binding energy. 

Another correlation that we were able to establish is 
between the maximum of the pair-binding energy and 
the position of the single-hole ground state in momentum 
space. In both the 4x4 and 6x6 periodic clusters the 
ground state shifts from the F — M and symmetry related 
directions of the Brillouin-zone to the zone-diagonals as t' 
is increased through tj^ax' ^^e Figs. [2]and[S] Specifically, 
exact diagnonalization^'* of the 4x4 cluster shows that 
the crystal momentum changes from (0,7r) and (tt, 0) to 
(0, 0) and (tt, tt) [CORE finds a similar transition to (tt, tt) 
but misses the (0, 0) state.] In the 6x6 cluster the shift is 
from (0,±27r/3) and (±27r/3,0) to (±27r/3, ±27r/3) [ex- 
cept for U = 1 — 3t where in a narrow region above ij„ax 




FIG. 10: The staggered magnetization in the two-hole and 
four-hole ground states of the 6x6 periodic chxster at U = 8f . 



the ground state is at (0, 0).] 

It is known from quantum Monte-Carlo simulations 
that the single-hole ground state of the homogeneous two- 
dimensional t— J model resides at (±7r/2, ±tt/2)}^ One 
may speculate whether this state is adiabatically con- 
nected to the ground state of the inhomogeneous model 
for t' > tjnax- The answer to this question is beyond 
the present study as it requires the diagonalization of 
larger clusters and the addition of higher-energy plaque- 
tte fermions with plaquette momentum (0,0) and (tt, tt) 
to the effective Hilbert space. Regardless of this point, it 
seems that the transition in the ground state momentum 
is a possible consequence of the maximum in Apb rather 
than its cause. We arrive at this conclusion based on the 
fact that in the 4x4 cluster Apf,(3/16) exhibits a max- 
imum of similar magnitude to that of Apb(l/16) while 
the 3-hole ground state is located at (0, 0) and (tt, tt) over 
the entire parameter rangeJ^ In the 6x6 cluster, on the 
other hand, the maximum in Apb(3/36) is accompanied 
by a change in the 3-hole ground state momentum, as 
depicted in Fig. [5l 



Phase stiffness 



In the thermodynamic limit of a d-wave superconduc- 
tor the pair-binding energy vanishes as Apt ~ 2AoA^~"'^/^, 
where Aq is the maximal value of the superconduct- 
ing gapj^ In our rather small clusters we can therefore 
roughly estimate Aq ~ Apb/ 2, which together with the 
d-wave BCS gap relation Tc — Ao/2.14, gives 



Tp = 



A 



pb 



(5) 



as a characteristic temperature at which pairs fall apart. 
The actual T^ may be smaller than Tp if phase fluc- 
tuations are important. To obtain an estimate for the 
phase-ordering temperature Tg we calculate the ground- 
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FIG. 11: The pairing scale T^ and the phase coherence scale 
Te in the two-hole and four-hole ground states of the 6x6 
periodic cluster at [/ = 8f. Te is shown for the cases where 
the phase twist is introduced at the bond level (solid lines) 
and at the plaquette level (dashed lines). The inset depicts 
Te of the two-hole system as the phase twist per bond (\) is 
varied from 7r/9 (upper curve) to 7r/72 (lower curve). 



state phase stiffness defines as 

1 d'^E 



Ps = 



A 



(6) 



</>=o 



Here E/A is the ground-state energy per unit area and 
(/) is a phase twist per bond in the x directioniiS. Ne- 
glecting the suppression of the stiffness due to thermal 
excitation of gapless nodal quasiparticles and using the 
relation Tc = 0.89ps for the two-dimensional XY model 
we obtain the estimator 



Te = Ps 



(7) 



We have calculated ps in two ways. In the first the 
phase twist was introduced into the Hamiltonian ^ by 
changing tij — ^ tije""^^-^ for two nearest-neighbor sites in 
the X direction. The effective CORE Hamiltonian for 
the twisted system was then derived and diagonalized 
to obtain the </> dependence of the ground state energy. 
In the second way the twist was introduced on the pla- 
quette level by modifying the couplings in the effective 
CORE Hamiltonian for the untwisted model ([T]). This 
was achieved via multiplication of a coupling between two 
neighboring plaquettes in the x direction that changes the 
number of holes on the right plaquette by An, by e*"^'^". 

The phase-ordering temperature of the periodic 6x6 
cluster with two and four holes is depicted in Fig. [11] 
The two methods yield similar results and they both en- 
counter problems in the region t' > t'^^^^. The nature 
of the difficulty is demonstrated by the inset of Fig. [TTl 
showing Ps as calculated from a discrete derivative of the 
ground state energy with respect to a twist introduced 
at the bond level. When the derivative is calculated for 



FIG. 12: The pairing scale Tp and the phase stiffness of a 14 x 2 
ladder in the two-hole and four-hole ground states at t/ = 
8t. Tp is essentially the same for the two hole doping levels. 
The solid (dashed) ps lines were calculated by introducing the 
phase twist at the bond (plaquette) level. 



increasingly smaller values of (j) the result does not con- 
verge for t' > t^jjax- Rather, it becomes negative and 
diverges, indicating that the CORE ground-state energy 
develops a cusp as function of (j). A similar behavior is 
also found in the 4x4 periodic cluster and in the lad- 
der systems. It occurs at lower values of t' for systems 
with odd number of holes. We take these findings as an 
indication that CORE is unable to produce a reliable ap- 
proximation for Ps in the region beyond the maximum in 
the pairing scale. 

In the range t' < t'^^^^ the estimated phase-ordering 
temperature increases monotonically with t' , but is con- 
sistently below the pairing scale. At t' = t'^;^^ we find for 
the two-hole system Tp/Tg ss 6. Increasing the doping to 
four holes decreases the maximal Tp slightly and increases 
Tg by about 70% leading to Tp/Tg{t' = C^^) « 3. The 
same holds true for the 4x4 cluster with two holes, which 
has a similar hole density and Tp/Tg{t' = ijjjax) ~ 2. Such 
a behavior suggests that superconductivity in the lightly 
doped two-dimensional checkerboard Hubbard model is 
governed by phase fluctuations. In ladders our definition 
Eq. ([6]) is equivalent to the phase stiffness along the lad- 
der (vcKc in the effective Luttinger liquid description of 
the system) divided by its width . As shown by Fig. [12] 
it is larger than the corresponding stiffness in the peri- 
odic clusters and grows with doping. However, since the 
one-dimensional system can not order it does not provide 
a phase ordering temperature similar to Eq. ([7]). 



D. Pairing correlations 

Another diagnostic tool for the presence of supercon- 
ductivity is the pair-field correlation function. We have 
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FIG. 13: Pair-field correlations at U — 8t in the ground state 
of (a) the 2-hole doped 4x4 periodic cluster, including a 
comparison to the exact results of Ref. 6,, and (b) the 2-hole 
and 4-hole doped 6x6 periodic cluster. D± and Dy are the 
correlations between the pair-field on bond ab and the pair 
field on bonds cd and ef, respectively, as defined by Fig. [T] 



calculated the following equal-time correlator 



ijM 



(Al,A«) 



(8) 



where ij denotes the bond between the nearest-neighbor 
sites i and j, and where the pair field on that bond is 
given by 



1 






itSi + St^ii 



"a)- 



(9) 



Fig. [13] shows the results for the pair-field correlations 
between the two most distant parallel (-Dii) and perpen- 
dicular (-Dj_) bonds on the periodic clusters with Nh, = 2 
and Nh — 4. Similar results were also obtained for the 
ladder systems. We find that Dy is positive and Dj_ is 
negative, consistent with d-wave pairing. The pairing 
correlations diminish in the limits t' /t — >■ and t' /t — ^ 1 
but unlike the pair-binding energy and the phase stiffness 
they are nearly independent of t' in the range of moderate 
inhomogeneity (from t' = O.lt to t' = 0.6i Apb,ps and D 
change by a factor of 7.5,4.5, and 1.5, respectively.) The 
magnitude of the correlations is small and comparable 
to results of previous studies of Hubbard ladders^ and 
Hubbard^ and t— J periodic clusters J^ 

The behavior of D suggests that pairing is very weak in 
the systems that were studied. This conclusion is in ap- 
parent contradiction with the large pair-binding energy 
found in the same clusters. In addition, as we already 



noted, the i'-dependence of the two quantities is very 
different. We believe that the fault may lie in the spe- 
cific form of the pair- field, Eq. ^, that was used for 
calculating the pairing correlations. It assumes a pair- 
wavefunction which is strongly localized in space. This 
may be wrong, as suggested by our results for the struc- 
ture of the ground-state. Fig. [HI clearly shows that most 
holes are not bound into pairs on a single plaquette. This 
is expected since for U — 8t the plaquette does not pro- 
vide a positive pair-binding energy. It seems, therefore, 
that thinking about Cooper-pairing in such systems in 
terms of real-space pairs occupying single bonds is a mis- 
leading oversimplification. Most likely, the phenomenon 
is more complicated and the pair wavefunction, while be- 
ing much more localized than its counterpart in conven- 
tional superconductors, still possesses a non-trivial real- 
space structure. 



IV. CONCLUSIONS 

This study had a dual motivation. First, to explore the 
utility of the CORE approximation as a method to in- 
vestigate fermionic strongly correlated systems, and sec- 
ondly to shed additional light on the role of inhomogene- 
ity in the physics of high-temperature superconductivity. 

As far as CORE is concerned, it is difficult to carry out 
the original scheme of Morningstar and Weinstein^ who 
iteratively applied the CORE method to obtain and an- 
alyze a fixed point Hamiltonian. In the case of the Hub- 
bard model there are simply too many couplings that are 
generated at each step. One is, therefore, forced to apply 
CORE once and investigate the resulting effective Hamil- 
tonian either by means of a mean-field approximatiork^, 
or via numerical diagonalization of finite clusters. The 
latter approach was previously implemented in the study 
of spin system o^^'^" and the t ~ J modeliSiii, and is the 
one which we pursued. As expected, when applied to the 
checkerboard Hubbard model range-2 CORE provides re- 
sults which are in good agreement with the available ex- 
act diagonalization results in the limit of small t' . In 
the moderate t' regime the method may be considered as 
semiquantitative and its validity in the uniform limit is 
questionable, particularly in the case of small U. More 
precisely, this statement depends on the property that 
one tries to calculate using the method. It seems that 
pairing is moderately local such that range-2 CORE is 
able to capture its salient features already in small sys- 
tems. The establishment of phase coherence, on the other 
hand, is a more extended phenomenon, for which the in- 
clusion of longer range effective interactions and diago- 
nalization of larger clusters are needed. In this context we 
would like to note that signatures associated with nodal 
quasiparticles of the putative d-wave Hubbard supercon- 
ductor, such as the suppression of the phase stiffness at 
low temperatures, are particularly difficult to capture us- 
ing range-2 CORE.^ 

Regarding the effects of inhomogeneity, our results 



10 



demonstrate that plaquettization of the Hubbard model 
may lead to a substantial enhancement of pairing. Opti- 
mal pairing is achieved at an intermediate scale of inho- 
mogeneity, which marks a crossover from a region with 
pronounced RVB characteristics to one with stronger lo- 
cal AFM correlations. The interactions of the doped 
holes with the spin background are the driving force of 
the pairing process. One should bare in mind, however, 
that the Hubbard plaquette, the building block of our 
model, is a special system. Its undoped ground state is a 
quintessential RVB state and it provides a positive pair- 
binding energy in a wide range of interaction strengths. 
Hence, it is interesting to ask whether a similar enhance- 
ment occurs for other plane patterns, especially those 
constructed from elementary clusters that do not exhibit 
pair binding. The possibility of such an outcome gains 
support from the fact that in the checkerboard model 
maximal pairing occurs at an interaction strength for 
which the pair-binding energy on each individual plaque- 
tte is negative. 

In the lightly doped clusters that we have studied su- 
perconductivity appears to be controlled by phase fluc- 
tuations. Owing to the reasons outlined above and our 
inability to carry out significant finite-size scaling it is 
difficult to estimate the phase ordering temperature in 
the two-dimensional limit and determine whether T^ in- 
deed achieves a maximum at an intermediate value of t' . 
Tc enhancement due to inhomogeneous pairing interac- 
tion was found in the attractive Hubbard model^i"— and 
the phase-ordering transition temperature is raised in the 
classical two-dimensional XY model with certain "frame- 
work" modulations of the phase couplings)^ We find it 
interesting to conclude by noting that Fig. [TT] hints at 
the possibility that a related inhomogeneity-induced en- 
hancement occurs in the model considered here as well. 
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Appendix A: The CORE Hamiltonian 

The full CORE Hamiltonian includes all possible terms 
that satisfy the symmetries of the problem, as detailed 
in Section |TT1 The resulting 45 effective couplings may 
be grouped in the following way 



The kinetic energy of the fermionic holes and the bosonic 
pairs is given by the first two terms. Ki^f contains 
the contribution of hopping processes involving only the 
charged degrees of freedom while K{,f+t contains simi- 
lar processes in which the triplet of AFM magnons also 
participate. The interactions among the fermions and 
pairs comprise VJ,/. Their remaining interactions with 
the magnon triplet, as well as couplings involving only 
the triplets, form the last group Vt- 

In the following, b], V^^ and fl^^^ create a hole pair, 
a magnon with spin component Sz = cr and a fermion 
with spin component Sz = <J and plaquette momentum 
q at site i, respectively. Our choice to use a basis where 
the two fermions have a definite plaquette momentum 
q — (0, tt) or q = (tt, 0) results in different interaction 
strengths between nearest neighbors in the x direction 
compared to the y direction. The notation {i,j)i, in the 
Hamiltonian below stands for nearest neighbors in the 
V = x,y direction and {AiBj)s^a signifies that the opera- 
tors Ai and Bj are coupled into an operator of total spin 
S and spin component Sz = cr. Finally, summation over 
S ,cr, q, and v indices is implied. 

The 7 "bare" kinetic couplings include fermion and 
pair hopping, as well as pair-fermion exchange and 
Andreev-like pair creation and disintegration. 






qo"2 



(ij). 



•^W E [6l/qtJqi.+4/qt./qi«+H.C.].(A2) 



(ij% 



Note that since the Hamiltonian is symmetric under ro- 
tations and reflections some of the couplings are related. 
For example, Jf ''^ — J?''', where q — q+ {tt, tt) mod 2tt. 
These symmetries ancl the d-wave symmetry of the pla- 
quette hole-pair state also imply J^;^ = —Jh'j. 

The remaining 9 kinetic couplings are associated with 
magnon-assisted hopping processes 



K, 



bf+t - Jbt > ,blt^ajb]iat 



+ j'r E itlfl,)sAt,u)s.o 



(^J% 



+ JfftT. [(44)i../q.. + H.c. 

(id). 



(A3) 



{^.3). 



H = Kbf + Kbf+t + Vtf + Vt. 



(Al) 



The 16 fermion and pair on-site energies and interac- 



tions are 

i i 

(id) {^,3). 



S,a 



(iJ) 



^2// ZJ ifcijfqi)s..<Tihtfqj)s.a 

(id). 



(A4) 



The fermion on-site energies depend on q in ladders 
where the symmetry between the x and y directions is 
broken. The on-site energies on a plaquette depend on 
the number of its nearest neighbors. Therefore, they may 
be position dependent in finite clusters without periodic 
boundary conditions. This does not happen for the clus- 
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ters that we have investigated. 

The last group consists of 13 couplings involving the 
magnons. They include their on-site energy, excitation 
amplitude from the vacuum, hopping matrix element and 
the strength of their mutual interaction together with 
their interaction couplings to the fermions and bosons. 
We find the coupling to the bosons to be very small. 



Vt - 



4^])o + H.c. 



(i-j) (i-j) 

^rE [(tl4)i,./q..+H.c.]. (A5) 
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